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Abstract 

Two integral relations derived from the Kohn Variational Principle (KVP) are used for describing 
scattering states. In usual applications the KVP requires the explicit form of the asymptotic 
behavior of the scattering wave function. This is not the case when the integral relations are 
applied since, due to their short range nature, the only condition for the scattering wave function 

is that it be the solution of (H — E)*$> = in the internal region. Several examples are analyzed 
for the computation of phase-shifts from bound state type wave functions or, in the case of the 
scattering of charged particles, it is possible to obtain phase-shifts using free asymptotic conditions. 
As a final example we discuss the use of the integral relations in the case of the Hyperspherical 
Adiabatic method. 
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I. INTRODUCTION 



The study of bound and scattering states in few-nucleon systems gives valuable infor- 
mation regarding the underlying nuclear interaction. The fact that the spectrum of each 
hydrogen and helium isotope has only one bound state in the mass region A = 2 — 4, limits 
the applicability of bound state methods to a few states in the study of these nuclei. In 
Ref. |l| a detailed study of the three-nucleon bound states has been done, whereas a similar 
analysis in the case of 4 He can be found in Refs. 

aa 

. In recent years much of the study 
of the three-nucleon system has been done in the three-nucleon continuum (see Refs. {4, ij] 
and references therein). Results in the four-body system have been obtained so far in the 
energy region below the three particle breakup p, |7|. 

Well established methods for treating both, bound and scattering states, regard the 
solution of the Faddeev equations (A = 3) or Faddeev-Yakubovsky equations (A = 4) in 
configuration or momentum space and the Hyperspherical Harmonic (HH) expansion in 
conjunction with the Kohn Variational Principle (KVP). These methods have proven to be 
of great accuracy. They have been tested using different benchmarks s[ Q . On the other 
hand, many other methods are presently used to describe bound states: for example the 
Green Function Montecarlo (GFMC) and No Core Shell Model (NCSM) methods have been 

nn 

used in nuclei up to A = 10 and A = 12 respectively [10J, [11[. Attempts to use these methods 



for the description of scattering states have recently appeared [12], Il3|. 

The possibility of employing bound state techniques for describing scattering states has 
always attracted particular attention [14]. Recently continuum-discretized states obtained 
from the stochastic variational method have also been used to study a + n scattering [l^. In 
those two approaches the tangent of the phase-shift results to be a quotient of two numbers. 
In the former the numerator and denominator are obtained from two integral relations after 
projecting the Schrodinger equation, whereas in the latter the numerator results from an 
integral relation derived by means of the Green's function formalism and the denominator 
from the normalization of the continuum-discretized state. 

Another problem that has received particular attention in few-nucleon scattering pro- 
cesses regards collisions between charged particles. Traditionally, the Faddeev method has 
been applied to the neutral n — d reaction. Applications to p — d zero energy scattering were 
studied in configuration space by the Los Alamos-Iowa group using s-wave potentials 
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and realistic forces [17| . In those calculations the KVP was used to correct the first order es- 
timate of the scattering length after solving the Faddeev equations in which the partial wave 
expansion of the Coulomb potential was truncated. Low energy p — d elastic scattering has 
also been studied using the pair correlated hyperspherical harmonic (PHH) expansion 18]. 
A benchmark comparing these two techniques was given in Ref. 19) . A different way to 
treat the Coulomb potential in few-nucleon scattering was proposed in Ref. 20], based on 



the works of Ref. 



2l|, in which the Alt-Grassberger-Sandhas equations were solved using a 



screened Coulomb potential and then the scattering amplitude was obtained after a renor- 
malization procedure. Summarizing, the description of scattering states using very accurate 
methods are at present limited to A < 4 systems. On the other hand, accurate methods 
for describing bound states beyond the A = 4 mass system exist, therefore the discussion of 
new methods for extending these approaches to treating scattering states is of interest. In 
this discussion the treatment of the Coulomb interaction cannot be neg lected. 

Recently two integral relations have been derived from the KVP 22J. It has been shown 
that starting from the KVP, the tangent of the phase-shift can be expressed as a quotient 
where the numerator and the denominator are given in term of two integral relations. This is 



similar to what was proposed in Ref. 



141 ]. however the variational character of the quotient 



and its strict relation to the KVP were not recognized. In fact, it is this property that makes 
possible many different and interesting applications of the integral relations. Accordingly, 
in the present study we would like to discuss some specific examples. We will show that the 
integral relations can be used to compute phase-shifts from bound state like functions. We 
start our analysis from the simplest case, the A = 2 system, using a model potential. Then, 
using a semirealistic interaction, n — d as well as p — d scattering are considered. This is 
of particular interest since, as we mentioned before, p — d scattering has been a subject of 
intense investigations. A second application of the integral relations regards the possibility 
of determining p — d phase-shifts from a calculation in which the Coulomb potential has been 
screened. Finally, as a third application, we will discuss the use of the integral relations with 
scattering wave functions obtained from the Hyperspherical Adiabatic (HA) expansion. All 
these examples serve to demonstrate the general validity of the KVP formulated in terms 
of integral relations. Due to their short-range nature, they are determined by the wave 
function in the interaction region and not from its explicit asymptotic behavior. This means 
that any wave function ^ satisfying (H — E)^ = in the interaction region can be used to 
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determine the corresponding scattering amplitude even when its asymptotic behavior is not 
the physical one. 

The discussion presented here is limited to systems with A = 2, 3. This is because 
our expertise to calculate few-nucleon wave functions is limited to these systems. However 
applications to heavier systems are possible and, in particular, it would be interesting to 
analyze the use of the GFMC method in the computation of the integral relations in systems 
with A > 4. The paper is organized as follows: in Section II the integral relations are derived 
from the KVP. Applications to the two- and three-body systems are given in Section III and 
IV, respectively, whereas applications of the integral relations in connection with the HA 
are given in Section V. The conclusions are given in the last section. 



II. INTEGRAL RELATIONS FROM THE KOHN VARIATIONAL PRINCIPLE 



In order to derive the integral relations we first consider a two-body system interacting 
through a short-range potential V(r) at the center of mass energy E in a relative angular 
momentum state 1 = 0. The solution of the Schrodinger equation in configuration space (m 
is twice the reduced mass), 

h 2 



(H - E)V(r) = ( V 2 + V - E)V(r) = , 

m 



(1) 



can be obtained after specifying the corresponding boundary conditions. For E > 0, with 
k 2 = E/(h 2 /m) and assuming a short-range potential V, ^(r) = 0(r)/v / 47r and 



0(r — > oo) — > Vk 
from which one gets \& — > AF + BG, where 

F = 
G = 



sm(kr) cos(A;r) 
A — ^ — - + B- v ; 




kr 



k sin(A;r) 



kr 



(2) 



4/T kr 
k cos(kr) 




(3) 



47r kr 

Use of the Wronskian theorem immediately leads to the following general expressions for 
the coefficients A and B: 



B 
A 



m 
m 



[< F\H - E\m > - < V\H - E\F >] 
[< m\H-E\G > - < G\H-E\V >] , 



(4) 



where we have made use of the fact that: 
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— [< F\H -E\G> - < G\H - E\F >] = 1 . (5) 

With the above normalization, and assuming that is an exact solution of Eq.([T]), it follows 
that ^ satisfies the following integral relations: 

m 

-— < V\H-E\F >= B 
n 2 
m 

— <^\H - E\G >= A (6) 
nr 

, B 
tano = — . 
A 

Explicitly they are 

Tfl i 

dr sin(kr)V(r)[r(p(r)} = B 

(7) 



o 



m f°° 0(0) 

dr cos(kr)V(r)[r(j)(r)] H -=■ = A 



H 2 \/k Jo Vk 
where in the last integral we have used the property that V 2 (l/r) = — A7i5(r). 

In practical cases the solution of the Schrodinger equation is obtained numerically. Then, 
tan 8 is extracted from 0(r) analyzing its behavior outside the range of the potential. The 
equivalence between the extracted value and that obtained from the integral relations defines 
the accuracy of the numerical computation. A relative difference of the order of 1CT 7 of the 
two values is usually achieved using standard numerical techniques to solve the differential 
equation and to compute the two one-dimensional integrals. The short range character of 
the integral relations should be noticed. This means that the phase-shift is determined by 
the internal structure of the wave function. 

The second relation of Eq. (J7J) shows a dependence on the value of the wave function at 
the origin. It could be convenient to eliminate this explicit dependence since the numerical 
determination of 0(0) might be problematic, as we will show. To this end we introduce a 
regularized function G = j reg G with the property that \G(r = 0)| < oo and G = G outside 
the interaction region. A possible choice is 

C=^f%-e- ), (8) 

where the regularization function f reg = (1 — e~ 7r ) has been introduced where 7 is a non 
linear parameter which will be discussed below. Values satisfying 7 > l/r , with r the 



range of the potential, could be appropriate. The regularized function G (as well as the 
irregular function G), satisfies the normalization condition 



m 



< F\H — E\G > — < G\H — E\F > =1. (9) 

Therefore the second integral relation in Eq ([6]) remains valid using G in place of G, 

in ~ 

— < m\H-E\G >=A , (10) 
nr 

with the explicit form: 



'in 



dr cos(kr)V(r)[r(f)(r)} + J 7 = A (11) 



h 2 Vk Jo 

where in J 7 all terms depending on 7, introduced by f reg , are included: 

1 f°° (m \ 
7 7 = — — — J dr y— V(r) cos kr — r ) 2 cos kr — 2^k sin krj e~ 7r [r0(r)] (12) 

Comparing Eq. f fTT]) to Eq. ([7]) we identify J 7 = (f)(0) j\fk. This equality can be verified with 
the same relative accuracy obtained for tan S provided that the regularization of G is done 
inside the interaction region. 

In the following we demonstrate that the relation tan 5 = B/A, which is an exact rela- 
tion when the exact wave function ^ is used in Eq (jH]), can be considered accurate up to 
second order when a trial wave function is used, as it has a strict connection with the Kohn 
variational principle. 

The connection of the integral relations with the KVP is straightforward. Defining a trial 
wave function fy t to be 

ty t = ^ c + AF + B G , (13) 

with ty c — >• as r — >■ 00, the condition \l> t AF + B G as r -> 00 is fulfilled. The KVP 
states that the second order estimate for tan S is 

[tan5] 2 " d =tan5-^ < (l/AW t \H - E\(l/A)% > . (14) 
n 2 

The above functional is stationary with respect to variations of \l/ c and tan 5. Without a 
loss of generality ^ c can be expanded in terms of a (square integrable) complete basis 

^c = ^a n (t) n (r) . (15) 



The variation of the functional with respect to the linear parameters a n and tan 5 leads to 
the following equations 

< <f> n \H-E\% >=0 

(16) 

<G\H-E\%>=0 . 
To obtain the last equation, the normalization relation of Eq. (Q has been used. From these 
two equations, \l/ c and the first order estimate of the phase shift (tan<5) is * can be determined. 
It should be noted that the first equation implies < ^f c \H — E\*ff t >— 0. Furthermore, from 
the general relation for A in Eq. @ , and using the second equation in Eq. (JTE]) , the following 
integral relation results 

77? ~ 

— <%\H-E\G>=A . (17) 
a 

Replacing the two relations of Eq.t flBT) into the functional of Eq. (fl4l) . a second order 
estimate of the phase shift is obtained 

[tan5] 2 " d = (tan5) ist - ^ < F\H - E\(l/A)% > . (18) 

Multiplying Eq. ( fl8l) by A one gets 

B 2 ' ld = B lst -^<F\H-E\%> . (19) 

On the other hand, a first order estimate for the coefficient B can be obtained from the 
general relation in Eq. (Tj0), i.e., 

— [< F\H - E\m t > - <%\H - E\F >]= B ls . (20) 

Therefore, replacing Eq. (l20l in Eq.( Tl9l . a second order integral relation for B is obtained. 
The above results can be summarized as follow 

3*"= ™ <% \ H -E\F> 
a 1 

A= ^ <^! t \H - E\G> (21) 
[tan5] 2nd = B 2nd /A. 

These equations extend the validity of the integral relations, given in Eq.([6]) for the 
exact wave functions, to trial wave functions. To be noticed that F, G are solutions of the 
Schrodinger equation in the asymptotic region, therefore (H — E)F — > and (H — E)G — > 
as the distance between the particles increases. As a consequence the decomposition of 
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ty t in the three terms of Eq. (fT3|) can be considered formal since, due to the short-range 
character of the integral relations, it is sufficient for the trial wave function to be a solution 
of (H — E)ty t = in the interaction region, without an explicit indication of its asymptotic 
behavior. This fact, together with the variational character of the relations, allows for a 
number of applications to be discussed in the next sections. 



III. USE OF THE INTEGRAL RELATIONS IN THE TWO-BODY CASE 

In this section we present applications of the integral relations of Eq. (|2ip to a two-body 



22 



23], we use a central, s-wave 



system. To make contact with the results given in Refs. 
gaussian potential 

V(r) = -V exp(-r 2 /r 2 ) , (22) 

with Vq = —51.5 MeV, ro = 1.6 fm and h 2 /m = 41.4696 MeV fm 2 . This potential has a 
shallow L = bound state with energy E-ib = —0.397743 MeV. 
We introduce the orthogonal basis 

TO = £g>(z)exp-(z/2) , (23) 

with C m a (normalized) Laguerre polynomial and z = 0r, where /3 is a nonlinear parameter, 
to expand the wave function of the system 

M-l 

*o=£ a m <p m . (24) 

m=0 

We solve the eigenvalue problem of H for different dimensions M of the basis. The variational 
principle states that 

E = (^o\H\^o) > E 2B , (25) 

with the equality valid when M — > oo. The nonlinear parameter (3 can be fixed to improve 
the convergence properties of the basis. In fact, for each value of M there is a value of 
f3 that minimizes the energy. Increasing M, the minimum of the energy becomes less and 
less dependent on j3 resulting in a plateau. Increasing further the dimension of the basis, 
the extension of the plateau increases as well, without any appreciable improvement in the 
eigenvalue, indicating that the convergence has been reached to a certain accuracy. At each 
step \l/o represents a first order estimate of the exact bound state wave function. 



S 



Since, in our example, the system has only one bound state, with appropriate values of 
M and (3, the diagonalization of H results in one negative eigenvalue Eq and M — 1 positive 
eigenvalues Ej (j = 1, M — 1). The corresponding wave functions 



are approximate solutions of (H — Ej)^j = in the interaction region. As r — > oo they go 
to zero exponentially and therefore they do not represent physical scattering states. The 
negative energy Eq and the first three positive energy eigenvalues (Ej, j = 1,3) are shown 
in Fig. [T] as a function of (3 for M = 40. We observe the plateau already reached by Eq for 
the values of ft showed in the figure. Furthermore, we observe the monotonic behavior of 
the positive eigenvalues towards zero as [3 decreases. The corresponding eigenvectors can be 
used to compute the integral relations of Eq. Q2ip and to calculate the second order estimate 
of the phase-shifts 5j at the specific energies Ej. This analysis is shown in Table H] in which 
the non linear parameter f3 of the Laguerre basis has been chosen to be 1.2 fm -1 . In the 
first row of the table the ground state energy is given for different values of the number 
M of Laguerre polynomials. The stability of Eq at the level of 1 keV is achieved already 
with M = 20. For a given value of M, Ej, with j = 1,2,3, are the first three positive 
eigenvalues. The eigenvectors corresponding to positive energies approximate the scattering 
states at these specific energies. Since the lowest scattering state appears at zero energy, 
none of the positive eigenvalues can reach this value for any finite values of M. We observe 
(see Fig. [T|) that the eigenvalues diminish as M increases. Defining fc| = p--£y, the second 
order estimate for the phase shift at each energy and at each value of M is obtained as 



On the other hand, as we are considering the A = 2 system, at each specific energy value 
Ej the phase shift tan<5j can be obtained by solving the Schrodinger equation numerically. 
The two values, [tan<y 2nd and tan<5j, are given in the Table [I] at the corresponding energies 
as a function of M. We observe that, as M increases, the relative difference between the 



Af-l 




(26) 




(27) 



[tan^f d = -^. 



i 



9 



variational estimate and the exact value reduces, for example at M = 40 it is about 10~ 6 . 
In fact, as M increases, each eigenvector gives a better representation of the exact wave 
function in the internal region and the second order estimates, [tan^] 2 ™ approach the exact 
result. 

The study of the stability of the results in terms of the non linear parameter 7 in 
the regularization function f reg is given in Fig. |2j In the upper panel, the second order 
[tan <5 1 ] 2nd , corresponding to the eigenvalue E\ given in Table HI is shown as a function of 
7, for M = 20,30,40. We observe a good stability for values of 7 > 0.2 fm -1 indicating 
that the regularization has to be done before w 5 fm. In lower panel the corresponding 
values for I 7 as a function of 7 are shown. The stable values obtained for M = 20,30,40 
when 7 > 0.2 fm -1 are I 7 = —8.6234, —8.4334, —8.3755 respectively. The exact values for 
4>(0)/y/k, obtained solving the Schrodinger equation numerically at the three energies are: 
—8.6188, —8.4338, —8.3755, respectively. We can observe that for M = 40 there is a com- 
plete agreement between I 7 and 0(O)/\/fc. Therefore J 7 can be considered to be an integral 
representation of (j)(0)/y/k. This is an important point since such a value can be used to 
normalize the variational wave function. In this example the integral relations derived from 
the KVP have been used to compute phase-shifts using bound state like wave functions. 

A different application of the integral relations regards the possibility of calculating the 
phase-shift of a process in which the two particles interact through a short range potential 
plus the Coulomb potential, imposing free asymptotic conditions to the wave function. As 
an example we use the same two body potential used in the previous analysis and add the 
Coulomb potential: 

e 2 

V(r) = -V exp -(r/r ) 2 + - . (28) 
For positive energies and / = 0, the wave function behaves asymptotically as 

¥ c) (r ->■ 00) = AF c (r) + BG c (r) , (29) 

with F c (r),G c (r) the regular and irregular Coulomb functions, respectively. The phase- 
shift is tan5 c — B/A. The KVP remains valid when the long range Coulomb potential is 
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considered and its form in terms of the integral relations results to be: 

-^<¥ t c) \H-E\F c >=B 



^<¥ t c) \H-E\G c >=A (30) 



<2?i d B 



[tan 5 ci - — 



with G c = f reg G c and a trial wave function behaving asymptotically as ty( c \ Since 
(H — E)\F C > and (H — E)\G C > go to zero outside the range of the short range potential, 
the integrals in Eq. (1301) are negligible outside that region. Therefore, for the computation 
of the phase-shift it is enough to require that \&j satisfies (H — E)^f^ = 0, inside that 
region. To exploit this fact, we introduce the following screened potential: 

^2 



V sc {r) = -V exp [-(r/r ) 2 ] + [e-^) n ] - . (31) 



For specific values of n and r sc it has the property of being extremely close to the potential 
V(r) of Eq. (1281) for r < ro, with r$ the range of the short range potential. The screening 
factor e~ (r / rsc ' )n cuts the Coulomb potential for r > r sc . Using the potential V sc to describe 
a scattering process, the wave function behaves asymptotically as 

^ n , rsc (r -> oo) = F{r) + tan 5 n , Tsc G(r) (32) 

where F, G are given by Eq. ( 1271) . since V sc is a short range potential. It should be noted 
that the screened phase-shift tan does not equal tan<5 c for any finite value of n and r sc . 
Solving the Schrodinger equation for V sc , it is possible to obtain the wave function ^ n ,r sc 
for different values of n and r sc . This wave function can be considered to be a trial wave 
function for the problem in which the Coulomb potential is unscreened. Accordingly it can 
be used as input in Eq. (1311 to obtain a second order estimate of the Coulomb phase-shift, 

< # nr \H-E\Fr >= B 



h 2 

TYl ~ 

-<m nr jH-E\G c >=A (33) 



r ei 2 nc ^ ^ 

[tand c J = — 

where the unscreened Coulomb potential is included in H. This estimate depends on n and 
r sc as the wave function does. In Fig. [3] the second order estimate [tan5 c ] 2 " d is shown as a 
function of r sc for different values of n. The straight line is the exact value of tan 5 C obtained 
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solving the Schrodinger equation. We can observe that for n > 4 and r sc > 30 fm the second 
order estimate coincides with the exact results. In this example the integral relations derived 
from the Kohn Variational Principle have been used to extract a phase-shift in the presence 
of the Coulomb potential using wave functions with free asymptotic conditions. 



IV. USE OF THE INTEGRAL RELATIONS IN THE THREE-BODY CASE 



The integral relations derived from the KVP are general and their validity is not limited 
to two-body systems. The two-body system is the simplest system in which different appli- 
cations can be studied and compared to the exact solution of the Schrodinger equation and, 
therefore, a detailed analysis of the variational character of the relations can be performed. 
In this section the study of the integral relations is extended to describe a 2 + 1 collision in 
the three-body system, below the breakup threshold into three particles. The description 
of the breakup channel remains outside the scope of the present work. In the following we 
will consider the two examples already discussed in the previous section: the computation 
of phase-shifts using bound state like wave functions and the calculation of phase-shifts in 
presence of the Coulomb potential using wave functions having free asymptotic conditions. 
To this end we will use the s-wave MT I-III nucleon-nucleon interaction 
singlet and triplet spin states, respectively: 



24j . active in the 



V M 1438 - 72 r - 3 .ii r 513.968 L55r 
Vmti{t) = e e 



V M- 1438 - 7 V 3 -^ 6 26.885 55r 

Vmtiii\?)— e e 

r r 



(34) 



with distances in fm and energies in MeV. This interaction has been used many times in the 
literature to study the three-nucleon system at low energies. It is considered a semi-realistic 
interaction since it describes reasonably well the deuteron binding energy and the singlet 
and triplet n — p scattering lengths. Its predictions for these quantities are Ed = —2.23069 
MeV, 1 a n _d = —23.582 fm, and 3 a n _^ = 5.5132 fm. To be noticed that this potential has a 
strong repulsion at short distances. 

To compute bound and scattering wave functions we make use of the pair hyperspherical 



harmonic (PHH) method which has proven to be extremely accurate 18|,|25). In the following 



a brief illustration of the method is given. For bound states, the three-nucleon wave function 
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is decomposed in three Faddeev-like amplitudes 

V = ip(X i ,Yl)+ip(X j ,Y i )+<tl;(X k ,Y k ) , (35) 

where we have introduced the Jacobi coordinates: X { = (rj — r k )/\/2 and Y, = (rj + 
r k — 2rj)/y / 6 (the generic vector indicates the position of nucleon k). Each amplitude 
having quantum numbers J,J Z ,T,T Z is expanded in angular-spin-isospin channels (called 
a-channels) as 

i/>{X u Y^ = ^ * a (Xi, Y^y a (jk, i) 

(36) 

and the radial amplitudes are expanded in terms of the PHH basis 

y<) = p'^/a^**) $>£(p) Wpk' 1 "*^) (37) 

K 

where we have introduced the hyperspherical variables, the hyperradius p and the hyperangle 
0i, defined by the relations X; L = p cos 4>i,Yi = p sin (pi, and ^ P^ ,La {4>i) is a hyperspherical 
polynomial. The summation is given in terms of the grand angular quantum number K = 
2n + l a + L a . The correlation functions f a ( r ) are introduced to accelerate the rate of 
convergence of the expansion. They take into account those configurations in which two 
nucleons are close to each other. A very convenient choice is to derive the correlation 
functions from a Schrodinger like equation governed by the two-body potential corresponding 



to the specific a-channel 25 |. 



In the following we consider a three- nucleon system in either the J = l/2 + or J = 3/2 + 
states with total isospin T = 1/2. The central character of the MT I-III interaction decouples 
those channels with different values of A a . Moreover, as the interaction acts only in the s- 
wave, we have l a = 0. This condition limits the number of channels of the (A Q = 0) J = l/2 + 
state to two channels, corresponding to = 0, 1, and to one channel, corresponding to 



s£ = 1, in the case of the (A a = 0) J = 3/2+ state. Finally, following Ref. |26j], the 
hyperradial functions are expanded in terms of Laguerre polynomials 

M-l 

u a K (p) = ^ m Ct\z) exp {-z/2) (38) 

m=0 

with z = (3p, and (3 a non linear parameter. We can define a complete antisymmetric three- 
nucleon state \a,K,m >, in terms of which the wave function \l/ n for the n-th state results 
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to be 

A a ^\a,K,m> . (39) 

a,K,m 

The linear coefficients in the expansion are determined by solving the generalized eigenvalue 
problem 

A i' n m' < a,K,m\H-E n \a',K',rri >= . (40) 

a',K',m' 

The extension of the PHH method to describe scattering states below the deuteron 



breakup, using the KVP, is straightforward 18|, |26|. As for bound states, we limit the 



discussion to the (A Q = 0) J = l/2 + ,3/2 + states with total isospin T = 1/2. The N — d 
scattering wave ^/ k function at the center of mass energy E = E d + (4/3)(h 2 /m)k 2 , is written 
as 

^k=J2 A K k m\ a > K > m> +I F * > +tan5 \G k > 

a,K,m 

\F k >=J2MX t )F (ky t )[s^l}jj z [P k ®l}TT z (41) 

i 

\G k >=J2MXi)freMGo(ky i )[s^^} JJ M k ® \W , 

i 

with 0rf(Xj) the deuteron wave function having spin s J = 1 and isospin P k = 0. F , Gq are 
proportional to the regular and irregular Bessel functions in the case of n — d scattering or to 
the regular and irregular Coulomb functions, divided by kyi, in the case of p — d scattering. 
The distance between the nucleon i and the deuteron, formed by nucleons j, k, is yi and 
freg(y) — (1 — ex P ( — iv)) is the chosen regularization factor. In our calculations value of 
7 = 0.25 fm -1 has been found to be appropriate. The coefficients A^ k m and the first order 
estimate of tan S are obtained by solving the following linear system 

A°^ k m , < a, K, m\H — E\at, K' , m! > + tan 5 < a, K, m\Q k >= — < a,K, m\T k > 

a',K',m' 

^ A% k m < a,K,m\g k > +tan5 < G k \Q k >= - < G k \T k > , 

a,K,m 

(42) 

where we have defined \Q k >= (H - E)\G k > and \F k >= (H — E)\F k >. Following Eq. ( 12TI) . 
the second order estimate for tan S is 



B 2 k = -j^<v k \T k > 



A k = ^<m k \g k > (43) 

[tan 4]^ = B 2 k nd /A k . 
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It should be observed that in the present case, due to the definition of the asymptotic 
behavior of ^k, we have (m/h 2 ) < ^k\Gk >= 1. 

In Table [Til t ne 3 H and 3 He bound states and the doublet and quartet n-d and p-d scat- 
tering lengths, corresponding to the MT I-III potential, are given in terms of the number M 
of Laguerre polynomials used in the expansion of the hyperradial functions. The calculations 
have been done using K = 16 which corresponds to 18 (9) hyperradial functions in the case 
of J = l/2 + (J = 3/2 + ). With M = 24, an accuracy better than 1 keV is obtained for the 
bound state energies and of the order of 0.001 fm for the scattering lengths. In Figs. H] and [5] 
the J — l/2 + ,3/2 + , I = 0, phase shifts 5 are given as a function of the energy in the form of 
the effective range functions for n — d and p — d, respectively. Following Ref. [27], for n — d 
scattering this function is defined as (E° = E — E^) 

K(E°) = k cot 5 (44) 

whereas for p — d scattering it is defined as 

K(E°) = C 2 { V )k cot 5 + 2krjh(r]) , (45) 

where rj is the Coulomb parameter, Cq = 27rr//(e 2?rT? — 1) and ^(77) = — ln(r/) + Reip(l + iff) 
(if) is the digamma function). The solid line in the figures represents these two functions 
obtained solving Eqs. (142 [ 43 1) for several values of the center of mass energy E° in the interval 
[0, \Ed\}. The solid points in the figures are the results obtained from the integral relations 
using bound state wave functions as explained below. 

The capability of the PHH to produce a very accurate description of bound and scattering 
states can be used to study different applications of the integral relations. The lowest 
eigenvalue after the diagonalization procedure of Eq. ( 140|) corresponds to the three-nucleon 
bound state energy of 3 H (T z = —1/2) or 3 He (T z = 1/2). However more negative eigenvalues 
could appear. For example, in the case of K = 16, M = 24 and (3 = 1 fm -1 , six negative 
eigenvalues E n appear satisfying |i? n | < \E^\. They are given in Table [TTTI transformed to the 
positive energies E® = E n — Ed- The corresponding eigenvectors approximately describe 
a scattering process at the center of mass energy E®, though asymptotically they go to 
zero. In the following we use the index n to label these approximate scattering states and 
reserve the continuous index k to label those scattering states having the correct asymptotic 
behavior, as given by Eq. fj4T|) . We now consider the diagonalization of the Hamiltonian 
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calculated using the PHH basis with the aforementioned values of K, M, and (3 but for the 
J = 3/2 + state. The J = 3/2 + state does not have any bound state, however six negative 
eigenvalues appear, all of them satisfying \E n \ < \E d \. The positive energies E® are also given 
in Table IIHI As in the previous case, the corresponding eigenvectors approximately describe 
the N — d scattering states, though asymptotically they go to zero. It should be observed 
that changing the values of K, M and 0, the number of these states and the corresponding 
energies at which the eigenvalues appear change. They do not present the stability that 
a true bound state shows. If we call E® the bound state energy calculated using a basis 
of dimension D, the variational principle establishes that Eq > E^b, with E%b the energy 
corresponding to D — > oo. When the value of D is sufficiently high, a further increase of the 
dimension of the basis will not give an appreciable improvement in Eq, showing a pattern 
of convergence of the type given in Table [TTJ On the other hand, the eigenvalues E n are 
embedded in the continuum spectrum of H which starts at E d . Accordingly, increasing D 
these eigenvalues tend to Ed and the number of them also increases. Similarly to what it has 
been done in the two-body case, we can consider these states to be approximate solutions 
of (H — E n )^ n = in the interaction region and use them as inputs in the integral relation 
to compute second order estimates of the phase-shifts. Defining k 2 n = (4/3)E®/(h 2 /m), with 
ty n the corresponding eigenvector and \F n >, \G n > the asymptotic functions of Eq. f HTj) 
calculated at k n , the second order estimate for the phase shift at each energy is obtained as 

--^ < * n \H-E n \F n >= B n 
2m 

— <y n \H-E n \G n >=A n (46) 
[tan<5 n r = ^. 

The second order estimates of the phase-shifts for the six cases given in Table IIHI are 
shown in Figs. H] and |5] as solid points in the effective range functions. We can observe 
an extremely good agreement with the scattering calculations. This method allows for the 
calculation of phase-shifts using bound state type functions, even in the case of charged 
particles. These results can be 

shifts were obtained solving the Faddeev equations in configuration space. For the n — d 



case, the results presented here and those from Ref. 2d\ are in complete agreement. In 



the p — d case the results of Ref. 271 were obtained considering the Coulomb potential in 

n 

s-wave, without including the correction obtained using the KVP as was done in Ref. |16| . 
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In fact, that paper reports the doublet and quartet p — d scattering lengths considering 
the Coulomb potential in s-wave (the given values are 0.16 fm and 13.75 fm, respectively). 



After the correction introduced 
potential, the results from Ref. 



using the KVP and considering the complete Coulomb 



16] are 0.003 fm and 13.95 fm, respectively. They are in 
close agreement with the results obtained here and given in Table [Til It is worth noting that 
the use of the integral relations permits a correct computation of the p — d phase-shifts in 
the energy range [0,Ey, after a diagonalization procedure of the Hamiltonian using square 
integrable basis functions. 

In the last example of this Section, we explore the possibility of extracting p — d phase- 
shifts from a calculation in which the Coulomb potential has been screened at a certain 
distance, as we have already done for the two body case. In the three nucleon system, we 
define the screened Coulomb potential as 

V sc (i,j) = [e-<'«/*--> n ] -(tl + l/2)(tj + 1/2) (47) 

with Tij the interparticle distance between nucleons Using the PHH method, we solve a 

p—d scattering problem using the screened potential and, therefore, the asymptotic behavior 
is of the form of Eq. (1411) . with Fq and Go the regular and irregular Bessel functions. For 
different values of n and r sc we calculate the scattering wave function \&™ ,Tsc and, using the 
integral relations, we determine the Coulomb phase shift. Similarly to what we have done 
in the two-body case, the integral relations are 

B l = -j^<K' 3C \H-E\F k > 

A k = ^ <^ rac \H - E\G k > (48) 

[tan5 c ,fc] 2 ' ld = Bf d /A k . 

In H the unscreened Coulomb potential is included and the asymptotic functions F k and 
Gk are given in Eq. (j4"Tl) in terms of the Coulomb functions Fq,Gq. The results are shown 
in Fig. [6] for the case of J = 3/2 + at E cm = 2 MeV. The second order estimates [tan5 Cj fc] 2 " d 
are given for different values of n as a function of r sc . For n > 4 and r sc pa 30 fm, the 
results are in complete agreement with the value obtained solving the p — d case without 
any screening of the Coulomb potential, tan<5 c = —2.1037, which is shown as a straight 



line in the figure. This method can be compared to the method given in Ref. 20) in which 



the Coulomb potential was screened using the same screening function as in Eq. fl47j) and 
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the Coulomb phase shift was recovered after a renormalization procedure. We can conclude 
that the integral relations used in Eq. f )48|) produce the same effect as the renormalization 
procedure. 



V. INTEGRAL RELATIONS WITHIN THE HYPERSPHERICAL ADIABATIC 
METHOD 



In order to study applications of the integral relations using the HA method for a three- 

3 (for 



nucleon system, we give a brief introduction to the method following Refs. 
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more details see Ref. 
expanded as 



28|). In the HA method the three-body wave function of Eq. (|35|) is 



(49) 



where $ M (p, fi) is a HA basis element and [p, ft] = [p, fa, Xi, Yj\ is the set of coordinates 
consisting of the hyperradius and of the five hyperspherical coordinates. The HA basis 
elements are the eigenf unctions of the hyperangular part of the Hamiltonian at fixed values 
of p: 



h 2 



2mp' 



: G 2 + v(p,n) 



%(p,n) = u IM (p)$,( P ,n) , (so) 

where G 2 is the grand-angular operator and V(p, Q) = J2i is the potential energy. The 

eigenvalues, t/ M (p), are the adiabatic potentials, that appear in the coupled set of differential 
equations 



&T P + U,(p)-£QM-E w,(p)- 



ri 2 



N A 



(51) 



with T p = d 2 /dp 2 + (5/ p)d/dp, the number of adiabatic channels included in the cal- 
culation, E the three-body energy, and from which the hyperradial functions w^p) are 
obtained. The coupling terms are defined as: P^i =< $ fi (p,Q)\d/dp\§ tl i(p,Q) > and 
Q^i =< <Iv(p, Vt)\d 2 /dp 2 \$pi(p, Vt)>. In Ref. Q the solution of the system of Eq. ([51]) has 
been studied for scattering states below the three-body breakup. In that study it emerged 
that the use of the integral relations helped to obtain a pattern of convergence for the 
phase-shift, in terms of the adiabatic channels, similar to the pattern obtained when the 
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HA expansion is used to describe the bound states. It was also shown that the convergence, 
without the application of the integral relations, is extremely slow. This problem originates 
in the boundary conditions imposed on the solution of the linear system. As p — > oo the 
scattering wave function behaves asymptotically as 



(f>d{r) 



sm(k p p) cos (k p p) 
tan dp — 



\ST). (52) 



with \ST > the total spin-isospin function and k 2 p = E°/(h 2 /2m) = \k 2 . In fact, when 
p — > oo, the distance Xj is limited by 4>d and yi = — > ^p, then the approximate relation 
ky-i ~ k p p holds. However, the exact equivalence between kyi and k p p is not matched for any 
finite value of p and, accordingly, the boundary condition of Eq. fl52|) is equivalent to that 
of Eq.(l4!l only at p w oo and Na — > oo. As a consequence 5 P converges extremely slowly 
to 5 by increasing the number of adiabatic states. Therefore, the application of the integral 
relations in the case of the HA method, as discussed in Ref. 22j, removes the limitation 
given by the slow convergence allowing an accurate description of the scattering states. 

Motivated by the results obtained in the previous section, we would like to analyze the 
possibility of computing phase-shifts solving the system of Eqs. (IST]) using bound state 
boundary conditions, namely imposing Wp(p) — > as p — > oo. To this end we expand the 
hyperradial functions in the basis of Laguerre polynomial given in Eq. ( 1381) and define the 
complete antisymmetric three- nucleon state \p, m >, with p indicating a HA basis element 
and with m a Laguerre basis element, respectively. In terms of this basis the three-body 
wave function results to be 

^n = Y, A U^ m > ( 53 ) 

where n indicates the different bound states. As we did in the previous section, fixing the 
number of adiabatic channels Na and the number M of Laguerre polynomials, we solve 
the generalized eigenvalue problem for specific values of J n ,T and identify the negative 
eigenvalues E n . They can represent true bound states (\E n \ > \E<i) or they can indicate the 
possibility of approximate scattering states (\E n \ < \Ed). As an example, for J = 3/2 + ,T = 
1/2, using the MT-III potential with AOi = 40, M = 60 and with the non-linear parameter 
fixed to the value (3 = 1 fm -1 , a very dense spectrum of 23 negative eigenvalues is obtained, 
all of them satisfying > \E d . The corresponding eigenvectors \l/ n are used to compute 
the second order estimates of the phase-shifts at the specific energies using Eqs. ( l4"6"j) . The 
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results are shown in Fig. [7] as solid points on the effective range line. We observe a perfect 
agreement between the 23 points and the exact results given by the straight line. In the 
solution of the same problem using the PHH basis, we have used M = 24 (see Table IIIip 
and obtained six negative eigenvalues. With the HA expansion, using M = 60, we obtain a 
much denser spectrum covering the whole energy range [0, Ey] below the breakup into three 
particles. 



VI. CONCLUSIONS 



The description of scattering states from the KVP has not been used in the literature as 
much as the equivalent form for bound states, the Rayleight-Ritz variational principle. A 
possible explanation for this is the different care that is required at the moment of describing 
the asymptotic structure of the system. For example, when a complete basis is used to 
describe an A-body bound state, the main condition for the basis elements is that they be 
square integrable. Elements having gaussian or exponential tails are often used. It is well 
known that these bases do not describe correctly the asymptotic structures as the distances 
between the particles increase. However the error introduced by these configurations in the 
computation of the binding energies is small. Conversely the extraction of the asymptotic 
constants could be problematic if the number of basis states is not sufficiently high (see for 



example Refs. 



3, 



29]). The situation drastically changes when the KVP is considered. The 



asymptotic structure of the system has to be introduced in an exact form in the trial wave 
function ty t otherwise the matrix element < *ff t \H — E\ty t > could well not be finite. In 
other words, \I/t has to satisfy asymptotically that (H — E)^> t = 0. Formally \I/t can be 
decomposed in an internal and in an asymptotic part as in Eq. (TIB"]) . Then, the internal 
part of can be expanded over a square integrable basis. However, the necessity of taking 
care explicitly of the (sometimes very complicated) asymptotic structures has limited the 
application of the KVP. The reformulation of the KVP, given in Eq. (|14p . to the form given 
in Eq. ([211) changes this situation. The KVP, given in terms of integral relations, does 
not necessitate the explicit introduction of the correct asymptotic behavior in \E^. The two 
integrals involved, < $? t \H — E\F > and < ty t \H — E\G >, go very fast to zero since F 
and G are asymptotic solutions of the Schrodinger equation. Therefore, the only condition 
necessary to obtain accurate second order estimates through the integral relations is that 
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the trial wave function fulfill (H — E)^t = in the interaction region. This condition can 
be achieved with a variety of methods. 

In the present paper we have discussed two applications of the integral relations: the 
use of bound state like wave functions to describe scattering states and, in the case of 
charged particles, the possibility of computing phase-shifts using scattering wave functions 
with free asymptotic conditions, obtained after screening the Coulomb interaction. Both 
problems are of interest in the study of light nuclei. We started discussing the applications 
to the A = 2 system with a model (short-range) potential. In this system the solution of 
the Schrodinger equation is possible and, therefore, meaningful comparisons between the 
variational estimates of the phase-shifts and the exact values can be performed. In the 
analysis it was shown that after a diagonalization procedure of the two-nucleon Hamiltonian 
those eigenvectors corresponding to eigenvalues embedded in the continuum spectrum can 
be used as inputs in the integral relations to determine the phase-shifts at those energies. 
We have observed that increasing the number of basis states, the phase-shifts converge to the 
exact values. In the second application we have performed a scattering calculation adding to 
the short-range potential a screened Coulomb potential. Accordingly we have imposed free 
asymptotic conditions to the wave function. It is well known that increasing the screening 
radius, the phase-shift calculated with the screened potential will never match the phase-shift 
obtained considering the full Coulomb potential. A renormalization procedure is necessary 
(see Ref. 2jJ and references therein). It is very interesting to observe that the relation 
integrals as given in Eq. ( 133]) produce the correct result. In fact, for suitable values of r sc 
and n, the wave function calculated with the screened potential, $fi™l, is an approximate 
solution of (H - E)vQ = in the region in which the short-range potential is active, with 
H containing the bare Coulomb interaction. In fact, due to the short-range character of 
the integral relations, it is equivalent to use *S or the wave function calculated with the 
Coulomb interaction in Eq. ( 133]) . 

These examples have been discussed also in the three-nucleon system. As a reference, we 
have used the PHH method which gives a very accurate description of the A = 3 system and 
is well documented in the literature. Firstly, we have calculated bound state wave functions 
using a semi-realistic interaction. For fixed values of J + and T the Hamiltonian has been 
diagonalized and attention has been given to those eigenvalues satisfying Ed < E < 0. This 
energy region corresponds to N — d elastic scattering and is located below the breakup into 
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three particles. The corresponding eigenvectors have been used to compute the second order 
estimate of the tangent of the phase-shift. In order to show the results in a visible way, we 
compute the effective range function K(E°) using the PHH method, which gives an almost 
exact result. Then, the second order estimates obtained from the bound state like wave 
functions have been compared to K(E°). We have observed that the variational estimates 
and the exact results at the level of four digits coincide. This is practically the level of 
accuracy reached by the PHH method, therefore we can conclude that the results based on 
the integral relations can reach the same level of accuracy as other methods usually used to 
describe scattering states in A = 3. Moreover, a similar accuracy has been obtained when 
the Coulomb interaction has been considered. We consider this result to be of particular 
importance. For example, the application of the Faddeev method for describing p — d 
scattering has been a subject of intense investigations and different techniques for including 
the long range Coulomb interaction has been proposed (see the Introduction). From the 
results presented here it emerges that elastic p — d scattering in the low energy region can be 
described using bound state like wave functions which can be easily computed. Furthermore, 
we have also analyzed the computation of p — d phase-shifts from a calculation in which the 
Coulomb potential has been screened, as we have done in the two-nucleon system. Again, for 
suitable values of r sc and n we were able to reproduce the p — d phase-shifts using the integral 
relations. This result will be useful for a simple extension of the Faddeev method, normally 
used to describe n — d scattering, to determine p — d phases without the normalization 
procedure described in Ref. j^. 

In the last Section we discuss an application of the integral relations in connection with 
the HA method. This method is often used to describe three-body bound states in nuclear, 
atomic and molecular physics. It is very efficient, in particular when the interaction has 
a strong repulsion at short distances. In Ref. 22] we have shown how to apply the HA 
method to describe a 1 + 2 collision solving the system of Eq. (I5TI) with appropriate boundary 
conditions and then, using the integral relations, the phase shift has been extracted. Here 
we have shown a different application, the system of equations given in Eq. ( |5T|) has been 
solved using bound state boundary conditions and then, using the integral relations, the 
phase shift has been extracted. We have solved the same case previously considered using 
the PHH expansion for the J = 3/2 + ,T = 1/2 state. We have obtained an extremely good 
description of the phase-shifts using HA bound state like wave functions. This application 
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will help to extend the applicability of the HA to describe, for example, atom-dimer collisions 
at low energies. 

Summarizing, we have demonstrated the usefulness of the KVP formulated in terms 
of integral relations. We have shown the general validity of the formulation with several 
applications to the A = 2, 3 systems. In particular, we have in mind the possible use of this 
technique for describing scattering states using bound state methods in systems with A > 4, 
as for example the GFMC method or the stochastic variational method. 
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M 



10 20 30 40 



E -0.395079 -0.397740 -0.397743 -0.397743 

Ex 0.536349 0.116356 0.048091 0.026008 

[tan^i] 2 "" -1.507280 -0.622242 -0.392005 -0.286479 
tan<5i -1.522377 -0.621938 -0.392021 -0.286480 
E 2 1.984580 0.449655 0.190019 0.103503 

[tan<5 2 ] 2nd -5.919685 -1.353736 -0.812313 -0.584389 
tan<5 2 -5.703495 -1.354691 -0.812270 -0.584388 
E 3 4.512635 0.994433 0.423117 0.231645 

[tan<5 3 ] 2nd 13.998124 -2.451174 -1.302799 -0.908128 
tan<5 3 12.684474 -2.448343 -1.302887 -0.908131 

TABLE I: The two-nucleon bound state Eq and the first three positive eigenvalues Ej (j = 1,3), 
as a function of the number of Laguerre polynomials M. The second order estimates, [tan<5j] 2nd , 
obtained applying the integral relations are given in each case and compared to the exact results, 
tan Sj . 
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M B( 3 H) B( 3 He) 


2 n 4„ 
""nd u nd 


2„ 
u pcf 


a-pd 


4 


-8.5117 


-7.8404 


1.0207 6.4590 


0.3987 


14.088 


8 


-8.5351 


-7.8683 


0.7251 6.4434 


0.0363 


13.978 


12 


-8.5357 


-7.8688 


0.7031 6.4413 


0.00636 


13.967 


16 


-8.5357 


-7.8689 


0.7019 6.4412 


0.00472 


13.966 


20 


-8.5357 


-7.8689 


0.7018 6.4412 


0.00461 


13.965 


24 


-8.5357 


-7.8689 


0.7018 6.4412 


0.00458 


13.965 


28 


-8.5357 


-7.8689 


0.7018 6.4412 


0.00456 


13.965 


32 


-8.5357 


-7.8689 


0.7018 6.4412 


0.00454 


13.965 



TABLE II: Convergence of the 3 H and 3 He bound states (in MeV) and the n — d and p — d doublet 
and quartet scattering lengths (in fm), using the PHH expansion, as a function of the number of 
Laguerre polynomials M. 
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3 H 



J = 1/2+ 
£°[MeV] [tan,5 n ] 2nd 


J = 3/2+ 
E°[MeV] [tan<y 2 " d 


0.05934 


0.03898 


0.06789 


0.30511 


0.18262 


0.09204 


0.20458 


0.54508 


0.39445 


0.17588 


0.43850 


0.84177 


0.70247 


0.28429 


0.77281 


1.21661 


1.11898 


0.41131 


1.21923 


1.74161 


1.65752 


0.55093 


1.79295 


2.59081 


3 He 


J = 1/2+ 
E°[MeV] [tan,5 n ] 2nd 


J = 3/2+ 
E°[MeV] [tan5 n ] 2 " d 


0.09857 


0.00753 


0.10338 


0.13807 


0.22822 


0.03915 


0.24470 


0.33970 


0.44642 


0.10965 


0.48326 


0.60099 


0.75903 


0.21248 


0.82144 


0.91555 


1.18003 


0.33871 


1.27157 


1.32466 


1.72398 


0.48111 


1.84927 


1.92215 



TABLE III: For the three-nucleon system, the six eigenvalues satisfying Ed < E n < (given in the 
form E® = E n — Ed) in the specific case of K = 16, M = 24 and (3 = 1 fm -1 . The corresponding 
second order estimate of [tan (5 n ] 2 " d , obtained from the integral relations, are also shown. 
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FIG. 1: (Color on line) The two-nucleon bound state energy Eq and the first three positive eigen- 
values Ej function of /3 in the case of M = 40 
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FIG. 2: The two-nucleon second order estimate, [tan5i] 2 " d , calculated using as a function of 
the non linear parameter 7, for the values M = 20, 30,40 and the integral l 1 as a function of 7 at 
the same three values of M. 



29 



-0.9 



CM 



CO° -0.95 
C 



tan5 



: 1 1 
\n=l 


l 

- 


\n=2 




n=3\ 


~~ 






n=5 n=6 . 

i i i i 


i 



10 



20 



sc 



30 

[fm] 



40 



50 



FIG. 3: (Color on line) The two-nucleon second order estimate [tan 5 c ] 2 " d as a function of r sc for 
different values of n. As a reference the exact value for tan<5 c is given as a straight line. 
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FIG. 4: The effective range function for J = l/2 + (a) and J = 3/2 + (b) in the n — d case. The 
solid points are obtained from the second order estimates of [tan5 n ] 2 " d given in Table lIlTl for 3 H, 
at the corresponding energies. 
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FIG. 5: The effective range function for J = l/2 + (a) and J = 3/2 + (b) in the p — d case. The 
solid points are obtained from the second order estimates of [tan(5 n ] 2n<1 given in Tabl dllll for 3 He, 
at the corresponding energies. 
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FIG. 6: (Color on line) The three-nucleon second order estimate [tan 5 C) k] 2 " d as a function of r sc 
for different values of n, at E = 0.2 MeV. As a reference the exact value for tan5 c is given as a 
straight line. 
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FIG. 7: The effective range function for J = 3/2 + (solid line) in the n — d case. The solid points 
are obtained from the second order estimates of [tan<5 n ] 2nd using the HA expansion 
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